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ABSTRACT 

Aims. We aim to study the evolution of dust attenuation in galaxies selected in the infrared (IR) in the redshift range in which they are 
known to dominate the star formation activity in the universe. The comparison with other measurements of dust attenuation in samples 
selected using different criteria will give us a global picture of the attenuation at work in star-forming galaxies and its evolution with 
redshift. 

Methods. We selected galaxies in the mid-IR from the deep survey of the North Ecliptic Field performed by the AKARI satellite. 
Using multiple filters of IRC instrument, we selected more than 4000 galaxies from their rest-frame emission at 8//m, from z ^ 0.2 to 
~ 2. We built spectral energy distributions from the rest-frame ultraviolet (UV) to the far-IR by adding ancillary data in the optical-near 
IR and from GALEX and Herschel surveys. We fit spectral energy distributions with the physically-motivated code CIGALE. We test 
different templates for active galactic nuclei (AGNs) and recipes for dust attenuation and estimate stellar masses, star formation rates, 
amount of dust attenuation, and AGN contribution to the total IR luminosity. We discuss the uncertainties affecting these estimates 
on a subsample of galaxies with spectroscopic redshifts. We also define a subsample of galaxies with an IR luminosity close to the 
characteristic IR luminosity at each redshift and study the evolution of dust attenuation of this selection representative of the bulk of 
the IR emission. 

Results. The AGN contribution to the total IR luminosity is found to be on average approximately 10%, with a slight increase with 
redshift. The determination of AGN contribution does not depend significantly on the assumed AGN templates except for galaxies 
detected in X-ray. The choice of attenuation law has a marginal impact on the determination of stellar masses and star formation 
rates. Dust attenuation in galaxies dominating the IR luminosity function is found to increase from z = 0 to z = 1 and to remain 
almost constant from z = 1 to z = 1.5. Conversely, when galaxies are selected at a fixed IR luminosity, their dust attenuation slightly 
decreases as redshift increases but with a large dispersion, confirming previous results obtained at lower redshift.The attenuation in 
our mid-IR selected sample is found ~ 2 mag higher than that found globally in the universe or in UV and Her line selections in the 
same redshift range. This difference is well explained by an increase of dust attenuation with the stellar mass, in global agreement with 
other recent studies. Starbursting galaxies do not systematically exhibit a high attenuation. We conclude that the galaxies selected in 
IR and dominating the star formation exhibit a higher attenuation than those measured on average in the universe because they are 
massive systems. Conversely UV selected galaxies exhibit a large range of stellar masses leading in a lower average attenuation than 
that found in an IR selection. 

Key words. Galaxies: photometry-Surveys-Infrared: galaxies-Ultraviolet: galaxies-Galaxies: evolution-Galaxies: ISM 


1. Introduction 

One of the most complex processes in galaxies is star formation 
and understanding it and its evolution with time remains a chal¬ 
lenge for modern astronomy. The evolution of dark matter is 
now modeled with high accuracy and is shaping the large scale 
structures in the universe. The physics of baryons, including star 
formation, is much more difficult to understand and to imple¬ 


ment in models of galaxy formation and evolution. The first nat¬ 
ural step is to find constrains from observations by measuring a 
reliable star formation rate (SER) at various scales and redshifts 
and studying the main drivers of its variations. 

Major progress has been made in the measure of the SER 
inside the disk of nearby galaxies with strong relation found 
between SER and molecular gas content (e.g., |Bigiel et al.|2008[ 
[Leroy et ar]|2008t [Kennicutt & Evans]|2012| for a review), and 


Article number, page 1 of 


15 









the measure of the star formation efficiency in galaxies is now 
possible up to large redshifts (e.g., Daddi et al. |201Q[ |Tacconi| 
|et al.pOlOl |2013t |Santini et al.| 2014| ). Besides these studies 
aimed to understand the process of star formation, the global 
amount of star formation in the universe is measured by building 
statistical samples with observables related to the recent star 
formation. The rest-frame UV emission is frequently privileged. 
The observations of the GALEX satellite up to z = 1 and the 
numerous optical surveys give very sensitive UV rest-frame 
measurements from the nearby universe to high redshifts and 
for large samples of galaxies (e.g., [Salim et al.|[2006t [Cucciati 


|et al.pOl^ [Finkelstein et al^|2012| . However , UV emission 
suffers from a major issue, which is dust attenuation. In the 
nearby universe and within galaxies like the Milky Way, about 
half of the stellar light is absorbed and re-emitted by dust at 
wavelengths larger than ~ 5yum, and that fraction can increase 
to more than 90% in starbursting objects. Therefore measuring 
the IR emission of galaxies has been identified as mandatory for 
measuring the total SFR, and adding both UV and IR emissions 
is now recognized as a very robust method to measure the SFR. 
The development of IR facilities allows the building of statistical 
samples of galaxies detected in IR, although the lower sensitivity 
of IR detectors combined with poor spatial resolution makes 
comparison with UV-optical surveys still difficult. The Spitzer 
and Herschel deep observations combined with ground-based 
optical surveys give us a complete view of the star formation up 
to z ~ 2 (e.g., [Reddy et al.|2008[ [Buat et al.|2012t [Reddy etaL| 

|2012[[oreo et al.|2014[) although selection biases are recognized 

to be complex ( [Bernhard et al.||2014t [Heinis et~ar]|2Q14| ). The 

very high redshifts (z » 2) remain almost unexplored in IR 
even with Herschel ^diddiU & Dickinson| ( 2014| ) and stacking 
methods are intensively used to find the average IR emission of 
optically faint objects not detected individually (e.g., [Reddy et 


al.|2012l|Ibar et al. |2013l|Heinis et al.|2013||2014l|Pannella et 

al.[2Q14[fwhitaker et al.|2Q14|r 

To complement the analysis of galaxy samples, the comparison 
of UV and IR luminosity functions gives the global energy 
budget of the universe. [Takeuchi et al.| ( |2Q05| ) first compare 
the UV and IR luminos ity functions from z = 0 to z = 1, and 
[Burgarella et al.| ( |2013| ) extend the analysis up to z = 4 taking 
advantage of the most recent Herschel measurements. Both 
studies measure the evolution with redshift of the global dust 
attenuation by comparing the UV and IR luminosity densities: 
they found it increases from z = 0 to z = 1.2 and then declining 
up to z=4 with a similar amount of dust attenuation at z = 0 and 
z = 4. Although these results depend on the exact shape of the 
IR luminosity function, which remains uncertain especially at 
high redshift ( [Casey et al.|2014a[ for a rev iew), they are found to 
be very consistent with those obtained by [Cucciati et al.| ( r2Q12| ) 
for a large sample of UV selected galaxies in the same redshift 
range. 

UV and IR luminosity functions exhibit very different shapes 
( [Burgarella et al.[[2013[ and references therein): bright objects 
are much more numerous in IR than in UV and the luminosity 
evolution with redshift is stronger in IR than in UV. The faint 
end of the IR luminosity function is not well constrained 
but may be considerabl y flatter than that mea sured for the 
UV luminosity function ( [Gruppioni et al.[[2013[ ). These very 
different distributions imply that UV or IR selected samples 
are expected to be composed of different galaxy population s, in 
terms of luminosity and dust attenuation ( [Buat et al.|2QQ7a[ ). 

The strong impact of dust attenuation and the need to estimate 
it as accurately as possible when UV data are used has led to 
numerous investigations of samples of star-forming galaxies. 
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Quite naturally these studies focus on galaxies selected either in 
their UV rest-frame or from their emission in gas recombination 
lines. A correlation is found between the amount of dust 
attenuation and the stellar mass of these galaxies at low and high 


redshift ([Garn & Best|2010 

[Sawicki 2012[[Buat et al. 

2012t[Ibar[ 

[et al. |2013l[Kashino et al. 

2013|l. Heinis etal.lpOM 

) study the 


IR properties of UV selected galaxies in the COSMOS field by 
stacking Herschel/SPIRE images at redshift 1.5, 3 and 4. They 
find that the average amount of dust attenuation correlates well 
with the stellar mass in their sample whereas it remains roughly 
constant when the observed UV luminosity (Luv) varies. For 
a given Lyv , the galaxies exhibit a very large range of masses 
and thus dust attenuation. 


The aim of the present work is to select galaxies in the 
IR to investigate the impact of this kind of selection on the 
dust attenuation properties and to compare our results to the 
ones obtained from UV rest-frame selections or globally in the 
universe. It is of particular importance since most of the star 
formation in the universe can be securely measured in IR up 
to redshift ~ 2 since the stellar emission reprocessed by dust 
dominates the unattenuated emission from stars. Numerous 
studies were already devoted to characterize dust attenuation in 
IR selected galaxies (e.g., Goldader et al. 2002}[Burgarella et"ar 


[2005l[Buat et al. [2005[[Howell et al. [2010|[Takeuchi et al.|2010[ ). 

A major result of these studies is the deviation of the sources 
detected in IR surveys from the relation found between the ratio 
of IR to UV luminosity and the slope of the UV continuum 
for local starburst galaxies with bluer colors than expected (but 
see also [Murphy et al.[ ( [2011[ )). Recently [Casey et ( [2014b[ ) 
attributed this departure to a patchy geometry. 

The infrared space telescope AKARI carried out a deep survey 
of the North Ecliptic Pole (hereafter NEP-deep) with all the 
filters of the InfraRedCamera (IRC). We take advantage of 
the continuous filter coverage in the mid-IR to build a 8yum 
rest-frame selection following the strategy of [Goto et al.[ ( [2010'| ). 
A 8 pm selection is relevant to select dusty galaxies active in 
star formation since it focuses on the emission of polycyclic 
aromatic hydrocarbons. Main sequence galaxies defined to have 
a normal mode of star formation and a tight relation between 
their star formation rate and their stellar mass are found to have 
similar ratios of Spm rest-frame to total IR luminosities ( [Elbaz[ 
[et al. [[20111 [Nordon et al.[[2012| [Murata et al.[[2014[ ). For these 
galaxies, a 8yum selection is similar to a selection based on their 
total IR emission. Using the Sll, L15, L18W, and L24 filters, 
we can select galaxies properly on a large range of redshift 
from z = 0.15 to z = 2.05. Combining the near to mid-IR 
cata log of [Murata et al.[ ( [2013[ ) with the optical data of [Qi et 
[ffi?[ ( [2Q14[ ) and with GALEX and HerscheljVACS detections, 
we are able to build the UV to IR spectral energy distributions 
(SED) of our selected sources. These SEDs are analyzed with 
the code CIGALE to measure physical parameters such as 
dust attenuation, SER, or stellar masses. We first optimize 
the SED fitting process on a subsample of galaxies with a full 
wavelength coverage and spectroscopic redshifts and then we 
run CIGALE on the whole sample of 4077 galaxies. The amount 
of dust attenuation and its evolution in redshift is measured and 
compared to the results found with other selections. 


The paper is organized as follows: the sample selection is 
described in Sect. 2, and the SED fitting process and its opti¬ 
mization on a subsample of galaxies is detailed in Sect. 3. The 
main result of this work, the measure of dust attenuation and 
its evolution with redshift, is presented in Sect. 4. In Sect. 5, 
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Table 1. Sample selection. Central wavelength and FWHM of the filters are expressed in jdm. The last column (UV data) is the number of galaxies 
for which at least one measurement in the UV continuum (<^ 0.25yum) is available. The filters corresponding to the UV range are indicated in 
brackets. 


Bin 

Eilters 

^cent 

EWHM 

Redshift range 

Number of galaxies 

PACS matches 

UV data 

1 

Sll 

10.61 

4.08 

0.15 <z<0.49 

1661 

352 

1031 (NUV) 

2 

L15 

15.98 

6.5 

0.75 <z< 1.34 

1906 

190 

1748 (u,g) 

3 

L18 

19.6 

11.2 

1.34 <z< 1.85 

460 

50 

409 (u,g,r) 

4 

L24 

23.11 

5.25 

1.85 <z < 2.05 

50 

7 

45 (g,r) 


we compare our results to other measures of dust attenuation for 
different sample selections and propose the stellar mass as the 
main driver for dust attenuation. Our conclusions are presented 
in Sect. 6. 

Throughout the paper we use the WMAP7 cosmological param¬ 
eters. The parameter Lir corresponds to the total luminosity 
emitted by dust as defined in Sect. 3, and the UV luminosity 
Luv is defined as vLy at 150 nm. All luminosities are expressed 
in solar units (L© = 3.8310^^ergs“^). 


2. Sample selection 

Our aim is to combine data at different wavelengths for galax¬ 
ies primarily selected at 8yum in their rest frame. We start with 
the Murata et al. ( 2013| ) catalog of the AKARI-NEP deep survey, 
which contains 27770 sources over 0.5 deg^. We cross-matched 
this catalog with the optical-near-IR catalog of |Oi et ar] ( |2014| ) 


for which photometric redshift are calculated. We found 23345 
sources in common within a tolerance radius of 1 arcsec on the 
coordinates. MegaCam (u*, g',r', i', z') and WIRCam (Y,J,Ks) 
data are available for these sources. The 5cr detection limits 
correspond to u^ = 24.6, g' = 26.5, r' = 25.7, i' = 24.9, z' = 
23.9, Y = 23.2, J = 22.8, and Ks = 22.5 mag (AB scale). The 
photometric red shifts are estimated using Le Phare software El 
bert et al. 12006]). The redshift accuracy (crAz/(i+ 2 )) is found to be 
equal t o 0.032 at z < 1 and 0.117 at z > 1. We refer to |Qi et alT 
( |2014| ) for a more detailed description. 

We perform our 8yum rest-frame selection following the same 
strategy as [Goto et al.| ( |2010] ). We consider the four AKARI fil¬ 
ters Sll, L15, L18 and L24 and define redshift bins so that the 
observed wavelength of the 8yum feature corresponds to a trans¬ 
mission of the filters larger than 0.8. We then check that the 
detection rate in the photometric catalog of |Qi et J^ ( [2014| ) is 
similar for each subsample of sources detected in one of these 
these filters (85% for sources detected in the Sll filter, 80% for 
the LI5 and LI8 filters, and 78% for the L24 filter), which en¬ 
sures us that we do not have any systematic bias between our 
bins. The LI8 filter has a very broad bandpass, which largely 
overlaps those of the L15 and L24 filters. We first define the red¬ 
shift range reachable with the L15 filter, and then the redshift bin 
of the LI8 filter starts at z = 1.34 in order not to overlap the one 
defined with the L15 filter. The L24 band has a low sensitivity 
and we prefer to use it only when the 8yum rest-frame comes out 
of the LI 8 filter (z > 1.85). We define a sample of 4077 galaxies. 
The redshift distribution is described in Table 1 together with the 
number of galaxies per bin of redshift defined with each filter. 
The AKARI-NEP field was also observed by GALEX for 89450 
s (Program GI4-057001-AKARI-NEP, P.I. M. Malkan). When 
available we add the GALEX/NGY data for galaxies with red¬ 
shift lower than 0.925, corresponding to a rest-frame wavelength 
larger than 1200 A. The NUV fluxes were measured from the 
images taken on the MAST database, the flux extraction is per¬ 
formed with DAOPHOT (Mazyed et al., in preparation) and the 


detected sources are cross-matched with the optical coordinates 
of the sources within a matching radius of 2 arcsec. NUV fluxes 
are added to the SED of 1207 sources with z < 0.925. 

The AKARI-NEP field was also observed by Herschel with the 
PACS and SPIRE instruments for 73.5 hours (Program OTl- 
sserjOl-1, P.I. S. Serjeant). Given the small number of sources 
detected with SPIRE in the field studied in this work, we only 
consider PACS data. Eluxes at 100 and 160 jam from the PACS 
images were also measured with DAOPHOT (Mazyed et al., in 
preparation) the match with our selected sources catalog is made 
with a matching radius of 3.5 arcsec, we keep only sources with 
a single match within this radius. We find that 599 sources are 
detected with PACS at 100 jam and 270, both at 100 and 160 jam 
(Table 1). 

We are interested in studying dust attenuation in these galaxies. 
To this aim, data in the UV range bring information very com¬ 
plementary to those in the IR range. Therefore we check the 
number of sources for which rest-frame UV data are available. 
In Table 1 we report the number of sources in each redshift bin 
with a least one detection at a wavelength lower than 0.25 [am in 
the rest frame of the galaxy. Approximately 90% of the sources 
are observed in UV except for the first redshift bin for which 
only 62% of the selected sources are detected with GALEX. 

The AKARI-NEP Deep Eield was observed with Chandra with a 
detection limit of ~ 10“^^ergs“^cm“^ in the 0.5-2 keV band. The 
X-ray source catalog is presented in |Krumpe et al.| ( [2015| ). We 
used the VLO version of the catalog. Unfortunately X-ray ob¬ 
servations cover only approximately half of our field. We cross- 
matched the catalog of X-ray sources from |Krumpe et al.| ( [2015| ) 
with our selection following their prescription (separation lower 
than the match radius) and found 147 matches. It corresponds to 
6% of our initial sample. We can expect the same low contribu¬ 
tion of such X-ray sources for the whole field. 

Last, we also define a subsample of our 8yum selected galaxies 
with a spectroscopic redshift from the compilation of |Qi et aT 


( |2014| ) and detected with PACS as well as in most of the other 
bands: 106 galaxies are thus selected, 7 are classified AGN type 
1, and 1 is classified AGN type 2. The field where spectroscopic 
redshifts are available is entirely covered by the X-ray observa¬ 
tions. Eifteen sources are detected, and four out of the seven 
AGN type 1 and AGN type 2 are also detected in X-ray. This 
catalog is called SPECZ sample in the following and is used to 
optimize the SED fitting process. The redshift distribution of the 
sources included in each sample is shown in Eigj^ 


3. Fitting the spectrai energy distributions 

The SED fitting is performed with the version vO.3 of the 
CIGALE code (Code Investigating GALaxy Emission J] devel¬ 
oped with PYTHON. CIGALE combines a UV-optical stellar 
SED with a dust component emitting in the IR and fully con¬ 
serves the energy balance between dust absorbed emission and 

^ http://cigale.lam.fr 
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Fig. 1. Redshift distribution of the whole sample (black solid line) 
of galaxies detected with PACS (black dotted line) and with a spectro¬ 
scopic redshift (SPECZ sample: green-filled histogram). 


Table 2. Range of values of the input parameters for dust and stellar 
emission, used for the SED fitting with CIGALE. 


Parameter 

Range 

Amount of dust attenuation E(B - Vy 

0.1-1 mag 

Attenuation curve 

B12,C00, SMC-like 

IR templates^, a 

1-3 

AGN fraction,fracAGN 

0-0.5 

Stellar populations 

Age (old stellar population) tf^ 

2-11 Gyr 

^-folding rate (old stellar population) r 

1-5 Gyr 

Age (young stellar population) tysp 

50-500 Myr 

Stellar mass fraction^ fysp 

0.01-0.2 


Notes. 

^ E{B - V) corresponds to the attenuation of the youngest population 
and a reduction factor of 0.5 is applied to the color excess of stellar 
populations older than 10^ years. 

^ Of is the exp onent of the powe r-law distribution of dust mass over heat¬ 
ing intensity ( [Dale et al.|2014| ). 

^ tf is always lower than the age of the universe at the redshift of the 
source. 

^ Stellar mass fraction produced with the young stellar population. 


its re-emission in the IR. This energy balance is also observed 
for the AGN emission absorbed and re-emitted by the dusty torus 
(see below). The total IR luminosity Ljr is defined as the sum of 
luminosities coming from stellar and AGN light re-processed by 
dust. Star formation history as well as dust attenuation charac¬ 
teristics, including the attenuation law, are input parameters that 
can be either taken free or fixed according to the available data 
or the specific aims. The main characteristics of the code are 
described in [Noll et~ST] p009| . We refer to Burgarella et al. (in 
preparation) and Boquien et al. (in preparation) for a detailed de¬ 
scription of the new version of the code. Here we only describe 
the assumptions and choices specific to the current study. The 
main parameters and range of input values are reported in Table 
2. The output values parameters are estimated by building the 
probability distribution function (PDF) and by taking the mean 
and standard deviation of the PDF. 
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3.1. Star formation history 

Different scenarios of star formation history are implemented in 
CIGALE: exponentially declining or rising SFR, delayed SFR, 
and a declining SFR with an overimposed burst. It is possible to 
study any SFR by adding a file. As shown in sever al studies (e.g., 
|Pforr et k.|2Q12[ Buat et al.|2014t [Conroy|2Q13[ and references 
therein), it is difficult to disentangle the different scenarios be¬ 
cause of the degeneracy of the fits. As long as IR data are avail¬ 
able and we are dealing with galaxies that are actively forming 
stars, the measure of dust attenuation is not affected by the de- 
tailed star formation history (jMeurer et al.|1999 [Buat et al.|2005t 


Cortese et al.|2Q08t|Buat et al.|2Qllt|Hao et al. We adopt 

a model with two stellar populations: a recent stellar population 
with a constant SFR on top of an older stellar population created 
with an exponentially declining SFR. Ages of the older stellar 
population and young component {tf and ^ysp respectively) are 
free parameters ( [Noll et al.|2009| ). The initial mass function of 
Chabrierj ( 2003) is adopted with the stellar synthesis models of 


Bruzual & Charlot|p003| ). Compared to the widely used, simple 


exponentially decreasing SFR, this type of bimodal star forma¬ 
tion history is able to better reproduce real systems with several 
phases of star formation known and to give more realistic stellar 
ages ( [Buat et al.|2014| and references therein). 


3.2. Dust attenuation recipe and dust re-emission 

Since we are fitting the SED of galaxies selected in the IR, 
the treatment of dust attenuation and re-emission is critical. In 
CIGALE, dust attenuation is controled by the choice of the at¬ 
tenuation curve. The light re-processed by dust and re-emitted 
in IR is modeled using IR SED templates. 


3.2.1. Attenuation curves 


To mo del the attenuation by dust, the code uses the attenuation 
law of jCalzetti et al.| ( [2000| ) (hereafter COO), and offers the pos¬ 
sibility of varying the steepnes s of this law and adding a bump 


centered at 2175 A. We refer to [Noll et al.| ( |2009| ) for a complete 
description of the dust attenuation prescription, which still holds 
in the new version of the code. In brief, the dust attenuation is 
described as 


A;t = E{B-V) {k'{A) + D;i^,y,E,{^)) 


Av 


( 1 ) 


where Ay = 5500 A, k'(A) comes from Calzetti et al. (2000 


Eq.4) and DxQ,y,E^{A), the Lorentzian-like Drude profile com¬ 
monly used to describe the UV bump ( [Fitzpatrick & Ma^ 
[1990[ ), is defined as 


2^,2 


D 




E^A^r 


( A ^ - ^ 2)2 + ^ 2^2 • 


( 2 ) 


The coarse wavelength coverage in the UV rest frame does not 
allow us to measure the dust attenuation curve for every galaxy. 
Instead we can run CIGALE with different predefined attenu¬ 
ation curves and compare the results. We obtained our fidu- 
cial attenuation curve by studying high redshift galaxies ( [Buat 
et al.[ ( [MT^ , hereafter B12) , corresponding to E’b = 1.6 and 


S = -0.27, which is close to the LMC2 extinction curve with a 
UV bump of moderate amplitude. We also consider the COO at¬ 
tenuation law (E’b = 0, b = 0). [Ilbert et ar] ( |2009[ ) found that con¬ 
sidering an SMC extinction curve and a screen geometry could 
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Wavelength(/inn) 

Fig. 2. SEP templates for AGNs considered in this work and defined 
from the [Fritz et ah] p006| ) library with different values of the optical 
depths T 9.7 and viewing angle if/. The other parameters are fixed as de¬ 
scribed in the text. Green line: T9.7 = 6 , = 0; red line: T9.7 = 1, = 0; 

blue line: T9.7 = 1, = 80; black line: T9.7 - 6 , 1 //- 80. The fiducial 

models correspond to T 9.7 = 1,6 and if/ - 0. The models represented 
here correspond to a bolometric luminosity of lO'^^ergs"^ 


help to measure photometric redshifts. The curve corresponding 
to b = -0.5 and £'b = 0 is similar to the SMC curve of |Gordon| 
|et al.| ( |2QQ3| ) ( [Boquien et al.|2012| ) and is added for comparison 
(hereafter SMC-like). We introduce a differential attenuation for 
young and old stars, the co lor excess is reduced by a factor 2 for 
stars older than 10^ years ( Calzetti||l997} Chariot & Fall||2Q00 
[Panuzzo et al. 12007] ). 


3.2.2. Dust emission templates 


The stellar emission absorbed by dust is re-emitted in IR. 
This emission is generated using different empirical templates. 
Draine & Li| ( |2007| ), |Casey| ( |?012| ), and |Dale et al.| ( |2014| ) models 
are implemented in CIGALE. The templates used in this work 
are those of [Dale et al.|(|2014| ). They are updated templates from 
the [Dale & Helou| ( ^02|) library, which provides good fits o f 
Herschel data (e. g., |Wuyts et al.||2Qllt [Magnelli et al.||2Q13] ). 
Dale et af] ( |2014| ) also add a quasar component. Here we use 
these templates without this quasar contribution since the AGN 
contribution is defined separately in CIGALE, as described be¬ 
low. 


3.3. AGN component 


With CIGALE, it is possible to add the contribution of an AGN 
to the SED. The adopted templates are those of |Eritz et al. 
( 2006| ) based on two components: the emission of the central 
source and the radiation from the dusty torus in the vicinity 
heated by the central source. [Eritz et af] ( |2006| ) introduce a set 
of six parameters in their model, describing the geometrical 
configuration of the torus and the properties of the dust and 
solving the radiation transfer equation. The dust distribution 
inside the torus is ass umed to be smooth. The diffe rent param¬ 
eters are described in |Hatziminaoglou et al.| ( |2008] ). When one 
accounts for this component in CIGALE, all the free parameters 
related to this AGN contribution are added to the already 


numerous parameters related to stellar (direct and re-processed) 
emission. With only photometric data we are confronted by the 
well-known degeneracy of SED fitting processes induced by 
the limited number of data points compared to the number of 
parameters used to build the SED templates. Erom simulated 


composite SEDs with an AGN and a stellar component, [Ciesla 
|et al.| ( |2015| ) find that only extreme values of the angle between 
the AGN axis and the line of sight, characterizing the nature of 
the AGN (type 1 or 2), can be estimated. Thus we choose to 
fix several p arameters to average values: from |Hatziminaoglou| 
|et al.| ( |2008| ) the ratio of the outer to inner radius is fixed to 
60, the dust density parameters p and y are fixed to -0.5 and 0, 
respectively, and the opening angle of the torus is equal to 100°. 
The remaining free parameters are the optical depth at 9.7 pm 
(T9.7 ) and the angle if/ between the AGN axis and the line of 
sight. 


CIGALE calculates the relative contribution of the dusty 
torus of the AGN to the total IR luminosity that is called the 
AGN fraction (fracACN)- We explore several combinations of the 
three parameters: T9.7 , ij/ , and fracACN- A large number of AGN 
fractions is introduced since it is the main parameter relevant to 
the AGN component that we want to retrieve: we consider eight 
values of fracACN linearly spaced between 0 and 0.5 (cf. Table 
2). Eollowing the prescription of [Ciesla et al.| ( [20T5] ), we con¬ 
sider only two extreme values for the angle = 0° and 80° to 
reproduce an AGN type 2 and type 1, respectively^. After sev¬ 
eral trade-offs we also consider a low 1 ) and a high ( t 9.7= 

6) optical depth model. [Hatziminaoglou et~ar] ( [2008| ) have also 
introduced models corresponding to low optical depth tori to im¬ 
prove the fits of their SWIRE/SDSS quasars. Our four selected 
templates are plotted in Eig|^ We can see that 6 and = 0° 
corresponds to a hidden AGN with a strong silicate absorption. 
The two templates corresponding io ij/ = 80° are very similar and 
correspond to a quasar emission with a strong UV/optical emis¬ 
sion. Combining this kind of an SED with a stellar component, 
even with a low AGN fraction, will induce a high contribution 
of the AGN. The template corresponding to 1 and = 0° 
exhibits a smooth distribution in the mid-IR, not very different 
from the cases with ij/ = 80° , if we ignore the small silicate 
emission of these templates, which does not affect broadband 
measurements. However, in this case, the UV-optical part is ob¬ 
scured since crossing the torus leads to a much lower flux than 
in the case if/ = 80°. The AGN models with a smooth dust dis¬ 
tribution are known to produce sudden jumps in the SEDs, and 
whether the line of sight intercepts the torus or not, a clumpy 
dust d istribution will produce a smoother transition ( [Eeltre et af 


2012[ ). So as not to introduce an overly large contribution of the 


AGN in the UV-optical, but to allow for a silicate absorption or a 
smoother distribution in the mid-IR, our fits will be based on two 
configurations: T9.7 = 1,6 and = 0 ° for scenario 1 , T9.7 = 6 
and = 0, 80 ° for scenario ^ 


^ The definition of ij/ used in the grid of AGN models implemented in 
the current v ersion of CIGALE (VO.3) is different from the one u sed in 
the paper of [Fritz et al.[ ( [2006[ ). If we note xj/Q the value used by [Fritz [ 
[et al.[ ( [2006] ) then if/ - 90 -xf/Q. For the sake of consistency, we use the 
CIGALE convention. 

^ We did not consider the four AGN templates in a same run to avoid 
putting too large weight ioxf/ - 80^, the two SEDs corresponding to this 
value of xf/ being very similar. 
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Fig. 3. Examples of fits for the SPECZ sample with various redshifts, values, and AGN contributions. The observed wavelength is plotted 
on the x-axis. On the y-axis the observed fiuxes are plotted with colored points. The best model and its different components (spectra of the 
unattenuated stellar, dust (thermal) and AGN emission) are plotted as lines. The different colors and lines are described in the legend inside the 
top left plot. In the top right plot the value of fracAON is not reported since we do not find any AGN component for this source. 


3.4. Fitting sources with spectroscopic redshifts 

The quality of the fitting process is checked on the SPECZ sam¬ 
ple for which spectroscopic redshifts are available. The fits are 
very good, 99 out of the 106 SEDs fitted with a reduced;^Hower 
than 3 and < xj >= 1-3. Some examples are given in Eigj^ 


3.4.1. AGN templates 


The SPECZ sample is well suited to check the influence of the 
choice of the AGN templates. We explore the two combinations 
of parameters presented above (scenario 1 and 2). In Eig. |^we 
compare the reduced and AGN fractions found with these two 
configurations, and X-ray sources are checked separately. In the 
majority of cases a good agreement is found for the;^^ and AGN 
fraction values. Only two values of T 9.7 or ij/ are introduced in 
both scenarios and the PDE of these parameters cannot be built: 
for the current discussion, we consider the values of T9.7 and if/ 
corresponding to the best model. When two input values of T9.7 
are considered (scenario 1, T9.7 = 1,6 ), the best model corre¬ 
sponds to T9.7 = 1 for 85 galaxies (80% of the cases). When two 
values of ij/ are allowed (scenario 2 , if/ = 0 , 80^), ij/ is found equal 
to 0° (AGN type 2) for the best model of 70 sources ( 66 % of the 
cases). Eight galaxies exhibit a difference between theirlarger 
than 0.5. All of these galaxies are detected in X-ray and three are 
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classified as AGN type 1 from their spectra. Their AGN fractions 
are also found to be more dispersed than for the other galaxies. 
In all but one case;^^ is significantly lower with scenario 1. Here¬ 
after we adopt scenario 1, corresponding to AGN type 2 models 
with either a high silicate absorption or an almost flat mid-IR 
continuum, as our fiducial scenario for the AGN component. 

It is worth noting that the choice of the AGN template substan¬ 
tially modifies its contribution to the UV-optical SED as shown 
in Eig|^ and we expect some impact on the determination of the 
attenuation. This impact is illustrated in Eig|^ with two exam¬ 
ples: a galaxy detected in X-ray and a galaxy defined as ’typical’ 
(without any X-ray emission or spectral classification). Their 
SEDs are fitted with the two scenarios inducing very different 
dust attenuations: 1.4 mag versus 6.4 mag for the X-ray source 
and 5.2 mag versus 6.8 mag for the ’typical’ galaxy. AGN frac¬ 
tions are similar and a low xl value is obtained for both scenar¬ 
ios. The two corresponding best models are shown: when only 
AGN type 2 templates are considered (left panels, scenario 1) the 
UV-optical SED is dominated by the stellar emission, whereas 
the quasar emission is found to dominate when the best model 
corresponds to an AGN type 1 (right panels, scenario 2). Ei- 
nally, we would like to emphasize the difficulty of fitting SEDs 
with multiple components. As shown in Eig|^ both fits are satis¬ 
factory with;^^ <1.5 whereas the nature of the AGN component 
found for each best model (type 1 or 2 ) differs. 








V. Buat et al.: Dust attenuation up to z - 2 in the AKARI North Ecliptic Pole Deep Field 




Fig. 5. Examples of fits with the two scenarios for the AGN component. Symbols and lines are the same as in Fig. Galaxy 1 (upper panels) 
is detected in X-ray and identified as an AGN type 1, galaxy 2 (lower panels) is a ’typical’ galaxy of our sample, without any X-ray detection or 
spectral classification. The left panels correspond to best models obtained with scenario 1 (T 9.7 - 1,6, xj/ - ) and right panels correspond to 

best models obtained with scenario 2 (T9.7 = 6 and if/ -0, 80®). All fits are good with consistent AGN fractions but the amount of dust attenuation 
varies according to the adopted scenario. 


3.4.2. Attenuation curves 


the V band and the lowest in the UV. 


Here we analyze the impact of changing the attenuation curve 
on the estimation of attenuation, SFR, and stellar masses. 
The code was run with the three different attenuation curves 
described in Sect. 3.2.1 (B12, COO and SMC-like). We adopt 
scenario 1 as our fiducial scenario for the AGN component 
(t 9 7 = 1,6 and 0^ = 0) and check that the results do not depend 
on the assumed AGN templates (scenario 1 or 2). The fits are 
very good for the three attenuations curves: in all cases more 
than 90% of the fits correspond to < 3. The SMC-like 
curve is chosen as the best choice for 25 sources against 38 
and 43 for COO and B12, respectively, implying a slightly lower 
pertinence for the SMC-like curve as compared to the two 
others. Systematic differences are found for Ajjy and Ay as 
shown in Fig|^ These differences occur because of the limited 
amount of data (if any) in the UV rest-frame range. Using the 
B12 curve as the reference, a larger attenuation is found in the 
V band for COO together with a lower attenuation in the UV, 
and inverse trends are found for the SMC-like curve. Opposite 
variations found for the UV and V bands are explained by the 
shape of the attenuation curve. The COO slope is flatter than the 
B12 and SMC-like slopes: for the same amount of IR emission, 
the attenuation is lower in UV and higher in the V band. The 
SMC-like curve is the steepest, giving the highest attenuation in 


Linear regressions between the different estimations yield 


Auv-coo - 0.89 X Auv-bi2 + 0.05, (3) 

Auv-smc = 1.12 X Auv-bi2 - 0-04, (4) 

Av-coo = 1.18 X Av-bi2 + 0.04, (5) 

Av-smc = 0-89 X Av-bi2 “ 0.03, (6) 


where the values obtained with the B12 curve (Av-bi 2 and 
AUV-B12) are taken as independent values for the regression. 

The SFR is mainly constrained by the IR emission, so we 
do not expect strong differences and only small variations are 
found. The variation of SFR is plotted against Auv in Fig|7] 
< A(log(SFR) >= -0.04 dex and -r0.02 dex between COO and 
B12, and SMC-like and B12, respectively. The discrepancy 
increases with attenuation as expected: the difference in attenua¬ 
tion increases with the attenuation itself. With the COO curve the 
UV attenuation is found to be lower, so is the SFR, the inverse 
holds for the SMC-like curve. 


The estimations of stellar masses are unaffected by chang¬ 
ing the attenuation curve as long as Ay <~ 1.4 mag. For higher 
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Fig. 4. Comparison ofxl and AGN fractions (fracAON) for the SPECZ 
sample. The results with scenario 1 for the AGN component (rgj -1,6 
and if / - 0”) are plotted on the x-axis and those for scenario 2 (T9.7 = 6 
and ij/ - 0, 80”) on the y-axis. Galaxies detected in X-ray, and galaxies 
spectroscopically classified as AGN type 1 are plotted with blue points 
and empty squares, respectively. The objects best fitted with either = 
or if/ - 80” for scenario 2 are differentiated (black points and pink 
empty triangles). The dispersion on the estimations of the AGN fraction 
(1 cr dispersion of the PDFs) is plotted with error bars. 


attenuations, when compared to the reference chosen here (B12 
law), the COO law gives larger stellar masses (corresponding to 
a larger dust correction in the V band) and the SMC-like curve 
lower masses (corresponding to a lower correction). The differ¬ 
ences remain lower than 0.1 dex on average (Fig|7]). 


4. Evolution of dust attenuation 

4.1. Definition of representative samples 


Our aim is to measure dust attenuation in galaxies selected with 
similar criteria at different redshifts. The combination of the dif¬ 
ferent IRC filters led us to perform a 8/im rest-frame selection. 
In Fig is reported the IR luminosity (obtained with CIGALE) 
against redshift for our sample. We estimated the detection limit 
at 5cr using the flux limits of |Murata et al.| ( |MT^ for each IRC 
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Fig. 6. Comparison of attenuation factors in the UV (upper panel) and 
in the V band (lower panel) for the different attenuation curves. The 
values found for the B12 attenuation curve are plotted on the x-axis, and 
those found for the COO (filled red circles) and SMC-like (blue empty 
circles) attenuation curves are plotted on the y-axis. The dispersion of 
the measures is reported as an error bar. 


filter used to define the sample (cf. Table 1), and the average ratio 
between Ljr (found by SED fitting) and the monochromatic lu¬ 
minosity in the IRC filter corresponding to the redshift selection. 
Clearly, we cannot study galaxies with a similar Ljr over the full 
redshift range since the limit in luminosity increases sharply with 
z. We adopt two strategies: 

1 . studying galaxies sampling the same domain of the luminos¬ 
ity function at different redshifts, 

2 . selecting galaxies of similar luminosity on a reduced range 
of redshift. 


In the first case, we must select the galaxies dominating the 
luminosity function and the luminosity density at a given red¬ 
shift. To this aim, we take the total IR luminosity functions of 
[Magnelli et al.| ( |2013|), who fitted a doubl e power-law function 
originally defined by [Sanders et al.| ( |2QQ3| ). In each redshift bin, 
we select galaxies with an IR luminosity inside a bin of 0.6 dex 
centered on the characteristic IR luminosity corresponding 
to the transition luminosity of the double power-law function. 
The lumino sity is calculated with relations given b }jMag- 
nelli et ah] ( [2013| ). The selection is represented in FigjSpana 


we split the bin 2 in two parts (z < 1 and z > 1). The selected 
sources lie above the 5cr detection limit for the first three bins, 
except for a small overlap at the end of the second redshift bin. 
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Fig. 7. Comparison of SFR and stellar mass (Mstar) estimations for 
different attenuations curves. The attenuation (in the UV (resp. V band) 
for the SFR (resp. Mstar) estimations) is plotted on the x-axis. The 
difference (expressed with a logarithmic scale) between SFR and Mstar 
estimations with the COO and the SMC-like attenuation curves and the 
B12 attenuation curve is plotted on the y-axis. The symbols are the 
same as in Fig§ 


The detection limit in bin 4 does not allow us to sample the 
luminosity function well, and we drop this bin for this selection. 
The average values of the stellar mass found for each redshift 
bin are re ported in Table. 3. They are close to the specific values 
found by [Ilbert et aL] ( [2010| ) for galaxies at similar redshift and 
with an intermediate star formation activity. These values are 
also in the range of values found by |Karim etaL] ( |2Q11| ) for the 
galaxies dominating the star formation rate density up to z = 3. 


13 
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Fig. 8. Redshift and luminosity distributions of the sources. The 
total IR luminosity Ljr, plotted on the y-axis, is an output of the SED 
fitting. The small dots represent the whole sample and the filled circles 
the subsample of galaxies selected around within each redsfhift bin 
(bin 1: red; bin 2: blue, split in two subbins corresponding to z < 1 and 
z > 1; bin 3: green; bin 4: magenta). The adopted values of log(Ljj^/Lo) 
10.84 (bin 1), 11.5 (bin 2, z < 1), 11.7 (bin 2, z > 1), 12 (bin 3), and 
12.35 (bin 4). The bin size in luminosity is 0.6 dex. Detection limits at 
5cr are reported as a black line. 


is lower than 5, with a similar result for each redshift bin and for 
galaxies detected with or without PACS. We are studying dust 
attenuation in the UV with a sampling of this wavelength range 
that changes with redshift. We must check that our results are 
not biased by varying wavelength coverage. We have also con¬ 
sidered the subsample of galaxies with available UV rest-frame 
data that we defined in Sect. 2 and Table 1. The results for this 
subsample are very similar to those found for the whole sample. 
The AGN fraction is on average approximately 10% for the 
whole sample and is found to globally increase with z (Fig|^. 
Our spectral coverage varies with redshift and, for most of the 
sources, we do not have PACS data. With only MIR data, the 
decomposition of the IR emission in a stellar and a nonthermal 
component is difficult. We checked that the same trend is found 
when only galaxies detected with PACS are considered as shown 
in Fig 1^ (as well as for our SPECZ sample, not shown here). The 
average value of frac agn is6+5,9+7,11 + 8,14+8, and 21 + 10% 
from bin 1 to bin 4. Scenario 2 (T9.7 = 6 and = 0, 80°) yields 
slightly lower AGN fractions but also increases with redshift (6, 
9, 11, and 15% from bin 1 to bin 4 with similar dispersions). Our 
results are in global agreement with those of |Goto et al.| ( |2Qll] ), 
who isolated AGN galaxies in their AKARI sample using 
emission line ratios and also found a very moderate contribution 
of AGNs to the global IR emission. 


In the second approach, we want to measure the evolution of 
galaxies of similar Ljr. We define a first luminosity bin, 11 < 
log(LiR/Lo) < 11.4, for galaxies with z < 1 corresponding to the 
bin 1 and a part of bin 2 (Figj^. At higher redshift, galaxies with 
Ljr > can be observed in the redshift bins 2 and 3, and 

we select galaxies in these bins with 11.7 < log(^iR/^o) < 12.1. 

4.2. SED fitting 

GIG ALE is run on the whole set of data (4077 sources) with a 
fiducial configuration (cf. Sect. 3). The reduceddistributions 
are very good: for 93 % of the sources the minimum value of;^^ 


4.3. Dust attenuation 

4.3.1. Evolution of the amount of dust attenuation with 
redshift for Ljj^ galaxies 

To study the average attenuation as a function of redshift, we first 
consider the subsample of galaxies with luminosities around L*^. 
The average values of Auv for this selection are reported in Ta¬ 
ble 3. 

We show in Sect.3 that the measure of the attenuation depends 
on the choice of the attenuation curve. We can calculate the aver¬ 
age values of the attenuation in each redshift bin for the different 
attenuation curves that we consider by applying regression rela- 
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redshift 

Fig. 9. AGN fractions measured with CIGALE plotted versus the 
redshift of the sources: the whole sample is plotted with grey points, 
the measures for galaxies detected with PACS are represented with red 
points. The average values and dispersion for each redshift bin are over 
plotted (black filled circles for the whole sample, red filled circles for 
the PACS detections). The AGN fractions are always larger than 0 since 
they are measured as the mean of the PDF. The last redshift bin is not 
considered for the PACS detections since it contains only 7 sources. 


tions (Eq. 3 and 4). The variation with redshift of the different 
estimations of the attenuation is shown in Fig[TQ| 

We must check that the definition of L*^ as well as the expected 
increase of the uncertainty on the measure of Auv with the AGN 
fraction do not affect our results. To this aim, we also consider 
the whole galaxy sample in each redshift bin and the sources 
with an AGN fraction lower than 10%. With the whole galaxy 
sample the average values of Auv are found 0.2 mag lower than 
for the subsample of galaxies with luminosities around L*^ ex¬ 
cept for the first bin for which the difference reaches 0.3 mag (as 
expected since we add a large number of fainter galaxies). When 
only galaxies with an AGN fraction lower than 10% are consid¬ 
ered, the difference never exceeds 0.2 mag. These variations are 
lower than the difference induced for example by changing of 
the attenuation law (Fig[TQ|. 


We can add other measures of attenuation, already obtained 
for IR selected galaxies, to our analysis. At z = 0, |Buat et| 
[3T| ( |2007a| ) combined GALEX and IRAS data and derived vol¬ 
ume averaged n ieasures of the atten uation. For a luminosity 


^iR ~ 10 ^^-^Lq ([Sanders et 
tion Auv = 2.4 + 0.7 mag 


1l| 


2003), the average UV attenua- 
uat et al.| p007b ) measured the 
attenuation of luminous IR Galaxies (LIRGs, Ljr > lO^^L©) at 
z = 0.7 by combining Spitzer and GALEX data and found a 
mean attenuation of 3.33 mag with a dispersion similar to that 
found at z = 0. |Choi et aL] ( |2006| ) also measured dust attenuation 
in mid-IR selected galaxies of the Spitzer First Fook Survey at 
z = 0.8 by comparing SFR measured with the strength of emis¬ 
sion lines and Ljr. We apply their relation between Av and Ljr 
to the average value of Ljr for the redshift bin 2 (Table 3) and 
get Av = 2.33 mag. The visual extinction can then be translated 
to an atten uation in the UV continuum as explained in |Buat et| 
[aL] ( |2007b| ) giving Auv = 3.4 mag. The dispersion is directly 
measured on Fig. 12 of |Choi et al.| ( 2006| ). Oi et al. (in prepara¬ 
tion) observe the Ha emission line of a sample of mid-IR bright 
sources at z = 0.9 with the Subaru/FMOS spectrograph. We 
consider their 25 confirmed detections: the average attenuation 
Ah^ = 1.5 + 0.60 mag, which translates to Auv = 3.4 + 1.3 mag 
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Table 3. Evolution of dust attenuation at UV wavelength with redshift 
for the fiducial model (B12 attenuation curve and AGN 2) and galaxies 
with a luminosity close to Mean values of redshift, Ljr, Mstar , 
and Auv as well as the number N of sources used are reported for each 
redshift bin (bin 2 is split in two subbins, see text). The dispersion of 
Auv is on the order of 1.3 mag in each bin (as plotted in Fig[T0|. 


Bin 

N 

< z > 

< log(LiR) > 
Lo 

log(Afstar) ^ 

Mo 

< Auv > 
mag 

1 

795 

0.4 

10.78 

10.41 

3.32 

2-1 

479 

0.9 

11.46 

10.61 

4.20 

2.2 

726 

1.2 

11.68 

10.73 

4.04 

3 

282 

1.5 

11.96 

10.96 

3.93 


(using the relations detailed in Sect.5.2). 

All these measures are overplotted on Fig[T^ Given the different 
methods applied to measure Auv they are found to be consistent 
and our obtained values extend the analysis up to z ~ 1.5. The 
attenuation in UV for Ljj^ galaxies is found to increase with red¬ 
shift from ~ 2.5 mag at z = 0 to ~ 4 mag at z - 1 and then to 
remain almost constant up to z = 1.5. 


4.3.2. Evolution of the amount of dust attenuation with 
redshift for galaxies of fixed Lir 

We also have access to the amount of dust attenuation in galaxies 
of similar Ljr on limited ranges of redshifts. For galaxies with 
11 < log(LiR/Lo) <11.4 and z < 1, we find < Auv >= 4 mag 
at < z >= 0.4 (bin 1) and < Auv >= 3.6 mag at < z >= 0.9 
(bin 2 with z < 1). The dispersion in both cases reaches 1.3 mag. 
For comparison, we can add the results of |Buat et~ar] ( |2007b| ) for 
LIRGs with < Auv >= 3.8 mag at z = 0 and 3.3 mag at z = 0.7 
with a dispersion about 0.7 mag. The results are consistent given 
the dispersion of the measurements and both show a slight de¬ 
crease of the attenuation when z increases. 

In the redshift bins 2 and 3, we select galaxies with 11.7 < 
log(LiR/Lo) < 12.1. The variation of Auv with z for these ob¬ 
jects is plotted n Fig[TT] Despite a very large dispersion, a de¬ 
crease of attenuation with z is found, the trend is also clearly 
found when only galaxies detected in their UV rest frame are 
considered. 


4.3.3. Variation of the amount of dust attenuation with Lir 
and Mstar 


Many previous studies have concluded that dust attenuation in- 
creases with SFR or bolometric emission of young stars (e.g., 


Hopkins et al.||2001l Martin et al.||2005l Choi et al. 2006[ Xu 

et al.||2007 Buat et al. 2007a[ Zheng et al. 2007) and stellar 

mass (e.g.,|Brinchmann et al.|2004 

Martin et al.|2007l Iglesias- 

Paramo et al. 

20071 IPannella et a 

l.||2009t IGarn & Best||2010l 

Buat et al.|20] 

[2| and reference therein). Recent studies find that 


the relation between atte nuation and stellar mass (Mstar) should 
not evolve with reds hift ( |Ibar et al~||2Q13| [Heinis et~ar]|2014t 
[Pannella et al.||2014| ). In Fig|12|' Auv is plotted versus Ljr and 
Mstar for our complete galaxy sample since there is no need here 
to restrict the analysis to any subsample. A general increase of 
Auv with Lir is found that explains the increase of the average 
attenuation with redshift for galaxies close to Ljj^. The very large 
dispersion leaves room for the slight decrease found with z for 
galaxies with a fixed Lir. 

When Auv is plotted against Mstar the redshift bins appear to 
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Fig. 10. Average values of the attenuation factor, Auv. versus redshift. The red symbols refer to IR selected samples. The results of the present 
study are plotted wi th red-filled circles fo r the fiducial modeling, empty circles and squar es are used for the COO and SMC-like attenuation laws. 
Measurements from |Buat et al.| ( [^07a|b| ) are plotted with red-filled triangles, those from |Choi et al.| ( [T006| ) an d Oi et al. with a filled red square 
and a filled red triangle. The blue symbols and lines refer to UV and Hq, lin e selections discusse d in Sect. 5. |Cucciati et al.| ( [20T^ results are 
reported with a blue-dashed line. The blue squares refer to measurements from Heinis et ^p014| ) based on a stacking analysis: the e mpty square 


Ibar et al. 


for the whole UV selection and the filled square for galaxies with Mstar = lO^^^M©. The blue circles refer to measurements from 
( 2013 1 for Hq, emitters: the empty circle is for th e whole selection and filled circle for galaxies with Mstar = IQ^Q ^Mq. The blue empty tr iangle 
refers to measurements from |Reddy et al.|(2008| ). The black line and crosses correspond to the global estimates of |Burgarella et al.| ( |2013| ) based 
on luminosity functions. Their analytical relation is plotted as a solid line and crosses correspond to thei r measurements at d ifferent redshifts. 
The dispersion of each measurement is reported with a vertical bar. The measure of the dispersion for the |Heinis et al.| ( [?014| ) values are based 
on a stacking analysis whereas the others are estimated b ased on individual detections. The green empty squares are the average values of the 
attenuation found with the model of [Bernhard et al.| ( f2014| ), discussed in Sect.5. 
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There is a global increase of Auv with Mstar with a very large 
dispersion. In FigfT2lwe overplot the relation between Auv and 

with a dispersion of 1 mag 
using Herschel detections. 
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Fig. 11. Attenuation factor, Auv. versus redshift for individual galaxies 
with 11.7 < log(LiR/Lo) < 12.1. Measurements for galaxies detected in 
their UV rest-frame are plotted with blue filled circles, the empty black 
circles correspond to galaxies without UV rest-frame data. 


overlap, except in the case of the first bin, which samples much 
lower masses than the others. No clear trend is found with z in 
agreement with a nonevolution with redshift at least for z > 0.5. 


Half of the sample, lying between the limits of the relation at 
+ 1(T and 35 %, are found above the upper limit, which roughly 
correspond to an excess of 20% of galaxies above the Auv^star 
relation of jHeinis et al.| ( |201^ and jPannella et al.| ( [2014| ). 


5. Discussion 

We now compare the redshift evolution of the attenuation 
measured with different methods and in samples selected 
with different criteria. We consider the global measure of 
[Burgarella et al.j ( |2013| ) based on luminosity densities as well 
as measurements performed in large samples of star-forming 
galaxies selected to be either UV or Ha emitters. We begin by a 
discussion of these measurements and then we compare them to 
our results for our mid-IR selection. 


5.1. Average dust attenuation in the universe and in UV 
selected samples 


[Burgarella et al.j ( [2Q13| ) measured the global attenuation in the 
universe by comparing the IR and UV luminosity densities. 
Their result is plotted in Fig[T^ (black line and crosses). 


Article number, page 11 of 15 































































bjO 

nj 

< 


0.15 < z < 0.49 
0.75 < z < 1.34 
1.34 < z < 1.85 
1.85 < z < 2.05 




. -i** 


lO 9.5 10.0 10.5 11.0 11.5 12.0 12.5 

log(L|R)(L0) 



IOg(M,t3r)(M0) 


Fig. 12. Attenuation factor, Auy, versus Ljr (upper panel) and M^tar 
(lower panel). The color coding is the same as in Fig[^ larger symbols 
are used for bins 3 and 4 to improve their visibility. In the lower panel, 
the relation of |Heinis et al.| ( [20T4] ) is plotted with a black solid line, the 
dotted line is the Icr dispersion of the relation (see text for details). 


This measure is not related to individual objects and can be 
considered as a measure of the attenuation averaged over all the 
galaxy populations. This measure is much lower than the one 
we measure in our sample selected in IR. The difference is on 
the order of 2 mag (1.6 mag for bin 1). 


Cucciati et al.| ( |2Q12| ) selected galaxies in their UV rest 
frame in the VVDS-0226-04 field up to z = 4.5. They derived 
dust attenuation for each galaxy of their UV selection through 
SED fitting (without IR data), the variation they found is in 
close agreement with the results of [Burgarella et al.] p013| ) as 
shown in Fig[T^ [Heinis et al.| ( |2013 ) measured the average 
attenuation of a UV selection at z ~ 1.5 in the COSMOS field 
by stacking Herschel/SFIRE images. The average value found 
for their whole selection is reported in Fig[T^ and is consistent 
with that derived by [Cucciati et al.| ( |2Q12|) at the same redshift. 
In Fig [T^ we also report the measure of [Reddy et al.| ( |2008| ) at 
z ~ 2 for a sample of galaxies selected on their UV rest-frame 
colors, which also agree with the other measurements. The 
amount of attenuation found for our IR selection is much higher 
than the values found either in a UV selection or globally in the 


Article number, page 12 of 15 



Fig. 13. Prediction of the average amou nt of attenuation in a UV 
selected sample at z=1.5 from the model of [Heinis et al.[ ( [20T^ . The 
threshold in Luv (observed) is plotted on the x-axis and the correspond¬ 
ing average amount of dust attenuation Auv is plotted on the y-axis with 
a blue line (see text for details). Th e dotted black horizontal line is the 
measure of [Burgarella et al.[p013[ ) at the same redshift. 


universe. 


At first glance the agreement between Burgarella et al. 
( 20I3[ ) and [Cucciati et al. \2012) or [Reddy et al. ( [2008[ ) mea- 
surements may be surprising since they are obtained using very 
different approaches. Cucciati et al. measured the attenuation in 
ind ividua l galaxies whereas the values reported by [Burgarella et 
aL] ( [2Q13[ ) and based on luminosity functions are averaged over 


the galaxy populations and cannot be directly compared to indi¬ 
vidual measurements. 

These similar results can be explained by the lack of correlation 
between the average attenuation and Luv observed, obtained by 
[Heinis et~ar] ( [2013| . We use their simulated catalog built to re¬ 
produce the IR properties of their UV selected sample at z = 1.5: 
we cut it at different values of Luv and calculate the luminosity 
density in UV and in IR of the selected objects. The ratio of these 
luminosity densities is then translated in Auv as in [Burgarella et 
[3T] ( [20131 for comparison purposes. The resulting attenuation is 
plotted against the threshold in Luy in Fig [13] and compared to 
the value found by [Burgarella et al.] ( [2QI3[ ) at the same redshift. 
It is clearly seen that the cut in Luv has no significant impact 
on the average measure of t he attenuation. Theref or e a galaxy 
sample selected in UV, as in [Cucciati et al.[p012[ ) or [Reddy 
[3T] ( [2008[ ), can be representative of the average attenuation in the 
universe, without systematic bias due to a limiting UV luminos¬ 
ity. 


5.2. Dust attenuation of Ha emitters 


[Ibar et af~[ ( [20I3[ ) analyzed the IR properties of galaxies de¬ 
tected in their Ha line (HIZEFS project) at z = 1.46. From a 
stacking analysis of Spitzer, Herschel, and AzTEC images they 
derived the average total IR luminosity of their sample. Using a 
local calibration between Ha line and IR continuum emissions 
from [Kennicutt et al. [ ( [2009[ ) they obtained a median attenua¬ 
tion Aua = 1.2 + 0.2 mag. This attenuation can be translated 
in UV using the recipe of [CalzettT[ ( [1997[ ). The attenuation Ahcf 


















































































































V. Buat et al.: Dust attenuation up to z - 2 in the AKARI North Ecliptic Pole Deep Field 



Fig. 14. Auv versus Mstar for starburst galaxies defined with a sSFR 
0.6 dex higher than the average value found in each redshift bin. The 
color coding is the same as in Fig[^ The Auv^star relation assumed 
by |Bemhard et al.| ( |2014| ) is plotted with a solid black line, and the Icr 
dispersion with dotted lines. 


is related to the color excess of the ionized gas by applying a 
standard extinction curve, i.e., 


EJB-y)=A^J2A5, 


(7) 


adopting a Galactic extinction curve ( |Pei|1992| ) and a screen ge¬ 
ometry. 

The color excess of the stellar continuum E{B-V) is related 
to the color excess of the gas ( |Calzetti|1997| ), i.e., 


E{B -V) = 0.44 X EAB - V). 


( 8 ) 


The attenuation in UV is then calculated with the adopted atten¬ 
uation law. To be consistent with our analysis, we use the B12 
attenuation curve. 


Auv = 12.56 X E{B -V) = 2.26 x Aho., (9) 

(Auv = 10.33 E{B -V) = 1.86 Ah^ if the COO law is used). 


The corresponding value, Auv = 2.7 + 0.4 mag (2.23+0.4 
for COO), is plotted in Fig[T^and is marginally consistent with 
the other measures performed in a UV selection or globally 
with the luminosity functions (and fully consistent if the COO 
attenuation law is used). [Kashino et al. | ( |2013| ) found a larger 
value for the ratio E{B - V)IEg{B - V) reaching 0.7-0.8. Using 
this value for the color excess would lead to an attenuation of 
4.6 mag, close to the values found for our IR selected samples, 
we do not favor this a high value as discussed below. 


5.3. Comparing dust attenuation in IR and UV or Ha selected 
samples 

We have to explain why the attenuation found in IR selections is 
much higher than those found in UV or emission lines selected 
samples and on average in the universe. With our IR selection, 
we select galaxies dominating the star formation at the redshifts 
considered. We find these galaxies experience a dust attenuation 
that is much higher than the average attenuation in the universe 
or that is found in UV selected samples: the galaxies dominating 
the star formation are not representative of the average attenua¬ 
tion in the universe. 


Dust attenuation is found to depend on the stellar mass without 
a strong evolution with redshift (see references above). We can 
try to use the relation linking these two quantities to account for 
the stellar mass distributions when we compare dust attenuation 
from different selections. 

Let us consider the measure of |Heinis et al.| ( |MT3| ) at z ~ 1.5: 
the average stellar mass for our IR selectio n in the redshift bin 


Heinis et al. 


(20141 


3 is lO^^-^M©. Using the relation found by 
between dust attenuation and Mstar, we find that the average at¬ 
tenuation increases from 1.8 to 3.8 mag if we only consider UV 
selected galaxies with an average stellar mass of lO^^-^M©, lead¬ 
ing to a value consistent with the measures obtained in IR sam¬ 
ples at the same redshift. 

The same exercise can be performed for the Ha selected sample 
of |Ibar et al~| ( |2013| ) at the same redshift. Using their relation be¬ 
tween attenuation in the Ha line and Mstar for Mstar = IO^^-^Mq, 
we find Alla = 1.6 mag, which translates to 3.6 mag in UV a 
result again consistent with what is found for our IR selected 
samples of similar mass. These new values, which account for 
the stellar mass distribution, are also reported i n FigfT0| Using 
the cal ibration of stellar and gas color excess of |Kashino et al. 
( |2013| ) gives an attenuation of 5.8 mag (and 4.7 mag for the COO 
law), much higher than that found in our IR selection. The value 
found by [Kashino et al. | ( [2013| ) may apply to galaxies experienc¬ 
ing a phase of intense star formation, but it seems to overestimate 
the UV attenuation in our sample of galaxies. However an inter¬ 
mediate value between 0.44 and 0.7-0.8 would also be consistent 
with our results. 

We find that considering the relation between dust attenuation 
and stellar mass helps at reconciling measures of dust attenua¬ 
tion in UV, Ha, and IR selections. The phenomenological model 
proposed by [Bernhard et al.j ( [2014[ ) to explain the evolution of 
UV and IR luminosity functions with redshift is a good frame¬ 
work to go ahead in our analysis. Starting from stellar mass 
distributions at different redshifts, they assume a SFR-Mstar re¬ 
lation (the so-called main sequence) and the relation bet ween 
attenuation and stellar mass found by [Heinis ctsl] p014[ ) with 
a significant dispersion (cr = 1 mag, see Fig[12[) to predict Ljr 
and Luv for each object. They are abl e to reproduce the a ver- 
age trends found in the UV selections of [Heinis et al.[ ( [2Q14] ) and 
their model is well suited to simulate an IR selection. We apply 
the same selection to their simulated sample as the selection we 
performed in Fig and Table 3. The average attenuation we get 
is equal respectively to 2.89, 3.93, 3.81, and 4.05 mag for the 

4 bins defined in Table 3, with a dispersion on the order of 1 
mag. It is in good agreement with our measurements as shown 

inFig[T^ 

A major assumption of the model of [Bernhard et al.[ ( [2014[ ), 
which we can try to check, is their recipe for dust attenuation. 
In their model they distinguish between galaxies on the main 
sequence and starbursting objects. The galaxies on the main se¬ 
quence follow the Auv “ Mstar relation described above and plot¬ 
ted in our Fig[T^ They are defined as starbursts sources with 
a specific SFR (sSFR defined as SFR/ Mstar) 0.6 dex above the 
average value of the main sequence. The attenuation of these 
starbursting objects is assumed to be higher than that of galaxies 
on the main sequence and the fraction of star formation attributed 
to the starburst is completely obscured in their model. We can 
check if starbursting galaxies from our sample exhibit a larger at¬ 
tenuation than the other galaxies. We adopt a similar definition 
to define starbursts in each redshift bin by selecting sources with 
a sSFR 0.6 dex above the average value of the sSFR found in the 
bin. Our aim is is not to fully study the relation between SFR and 
Mstar in our sample, and moreover our S/tm selection is likely to 
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miss starbursting systems with an IR luminosity close to the de¬ 
tection limit since these objects exhibit a lower 8//m rest frame 
to total IR luminosity ( [Murata et al.|2014| ). The stellar mass and 
dust attenuation of these galaxies with a large sSFR are reported 
in Fig[^ We find that the starburst galaxies are preferentially 
located above the average relation, thus exhibiting a large at¬ 
tenuation, but the effect is not systematic. These galaxies are 
neither numerous nor attenuated enough to explain the ’excess’ 
of ~ 20% of galaxies found in Sect. 4.3.3 to show an attenuation 
larger than that expected from the average Auv - Afstar relation: 
in our sample, galaxies with a more moderate star formation ac¬ 
tivity can also experience a large dust attenuation. 


6. Conclusions 

We have studied the evolution of the amount of dust attenuation 
in a sample of 4077 galaxies selected in their rest-frame 8/im 
in the AKARI-NEP Deep Field observed with the IRC instru¬ 
ment. To complement optical and near-infrared data, GALEX 
and Herschel/FACS data are also added to the set of data. 

Dust attenuation is measured from the fit of the UV-to-IR SED 
of each galaxy with the CIGALE code. An AGN component is 
introduced from the Fritz et aT] ( |2006| ) library. We use differ¬ 
ent attenuation laws. The method of SED fitting is optimized 
on a subsample of 106 galaxies with spectroscopic redshift. The 
choice of the attenuation law is found to affect only marginally 
the measure of stellar mass and of SFR, the latter being essen¬ 
tially measured with the IR emission. Conversely, the amount of 
dust attenuation in the UV rest frame at 150 n ms varies with the 
choice of the attenuation curve: adopting the |Buat et al.| ( |2012 ) 
law as a reference, a higher attenuation is foun d with a SMC 
extinction law and a lower attenuation with the ICalzetti et al.l 


(2000) recipe, and inverse trends are found for the attenuation in 
the rest-frame V-band. We find that the AGN contribution to the 
total IR emission is low, but the introduction of AGN type 1 or 
type 2 templates in the fitting process affects the estimation of 
dust attenuation. The dust attenuation is systematically higher 
with an AGN type 1 contribution. 

Dust attenuation is found to increase with redshift up to z = 1 
and the measurements are consistent with a constant attenuation 
in the redshift range 1-1.5, in agreement with previous measure¬ 
ments on a smaller range of redshift. Dust attenuation increases 
with both Lir and Mstar- No trend with redshift is found when 
Auv is plotted against Mstar- For a given Ljr, we find the dust 
attenuation decreases slightly when redshift increases in consis¬ 
tency with the large dispersion found when the attenuation is 
plotted against Ljr. We find that the dust attenuation is larger 
than that measured in UV selected samples or samples selected 
in the Ha line at similar redshift, and that the dust attenuation is 
larger than the global attenuation in the universe measured with 
the IR and UV luminosity densities, the difference reaching ~ 2 
mag. This large difference is well explained when considering 
the average relation between dust attenuation and stellar mass: 
an IR selection is biased toward massive galaxies whereas UV 
selected galaxies exhibit a large range of stellar mass. Starburst¬ 
ing galaxies do not systematically exhibit a large attenuation. 
The phenomenological model of [Bernhard et al.| ( [2014| ) is able 
to very well reproduce the average measure of dust attenuation 
at different redshifts, but we find that the relation between Mstar 
and Auv more dispersed in our sample than is assumed in their 
model. 
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